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ABSTRACT 

We present a method for computing uncertainties in spectral models, i.e. level 
populations, line emissivities, and emission line ratios, based upon the propaga- 
tion of uncertainties originating from atomic data. We provide analytic expres- 
sions, in the form of linear sets of algebraic equations, for the coupled uncer- 
tainties among all levels. These equations can be solved efficiently for any set 
of physical conditions and uncertainties in the atomic data. We illustrate our 
method applied to spectral models of O III and Fe II and discuss the impact of 
the uncertainties on atomic systems under different physical conditions. As to 
intrinsic uncertainties in theoretical atomic data, we propose that these uncer- 
tainties can be estimated from the dispersion in the results from various indepen- 
dent calculations. This technique provides excellent results for the uncertainties 
in A- values of forbidden transitions in [Fe II]. 

1. Introduction 

Much can be learned about the behavior and evolution of an astronomical source 
through analysis of properly calibrated spectra. Interpretation can lead to estimates of 
density and temperature conditions, the chemical composition, the dynamics, and the 
sources of energy that power the emitting object. Such an interpretation requires modeling, 
with sufficiently high accuracy, the excitation and ionization balance of plasmas out of local 
thermodynamic equilibrium (LTE). But ultimately, the accuracy of the models depends on 
the quality of the atomic/molecular data employed. 

At present, atomic data exists for most spectral lines observed from the infrared 
to the X-rays. These data account for most processes leading to tens of thousands of 
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transitions from all ionic stages of nearly all elements of the first five rows of the periodic 
table. However, this huge amount of data has been obtained primarily through theoretical 
calculations with only sparse checks against experimental measurements. 

Despite many advances in spectral modeling, mostly in terms of increased completeness 
and improved quality of atomic/molecular parameters, a generalized quantitative estimate 
of uncertainties of the resultant models has not been provided. In recent years, a few 
authors have presented methods based on the Monte C arlo numerical technique for 



propa g ating uncertainti e s thro ugh spectral models, e.g. 



(|2Q12|); 



Ballance et al. 



Wesson. Stock, and Scicluna 



(120131 ) . As these techniques are very inefficient, their general 
applicability to complex spectral modeling is very limited. Finding a general and efficient 
method for estimating uncertainties in spectral models is important for two reasons. 
First, the accuracy of atomic/molecular data must be known before reliable conclusions 
can be provided on physically realistic comparisons between theoretical and observed 
spectra. At present researchers can only provide best fits to observed spectra without 
much understanding of the uncertainties impacting the results. Second, homogeneously 
accurate atomic data for all transitions of a complex and/or very large atomic system, like 
for example systems with multiple metastable levels (e.g. Fe II) or models with hundreds 
of energy levels as those needed in UV and X-ray spectroscopy, cannot be obtained. In 
such models, error propagation analysis of the spectrum could discriminate between a few 
critically important atomic transitions and the very large numbers of less consequential 
transitions. Conversely, detailed error analysis could direct further theoretical and/or 
experimental efforts to selectivelly obtain specific atomic measures that would significantly 
improve spectral models, rather than trying to determine all possible rates at once. 

This paper is organized as follows. Section 2 presents the analytical solution to the 
uncertainties in level populations of a non-LTE spectral model for assumed uncertainties in 
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atomic parameters. In Section 3 we propose a mechanism to estimate the uncertainties in 
atomic/molecular data and we test this by the case of Fe II through extensive comparisons 
with observed spectra. In Section 4 we discuss the uncertainties in line emissivities and 
emission line ratio diagnostics. Section 5 presents our conclusions. For the sake of clarity, 
the rest of the paper deals explicitly with the case of population balance by electron impact 
excitation followed by spontaneous radiative decay. It is also assumed that the plasma is 
optically thin. However, we note that our method can easily be extended to ionization 
balance computations, to additional excitation mechanisms such as continuum and Bowen 
fluorescence, and to optically thick transitions. 



2. Uncertainties in Level Populations and Column Densities 

Under steady-state balance the population, N i: of a level % is given by 

where n e is the electron density, A^i is the Einstein spontaneous radiative rate from 
level k to level % and q^j is the electron impact transition rate coefficient for transitions 
from level k to level i. Here, we assume that the electron velocity distribution follows the 
Maxwell-Boltzmann function, thus q^^ and q^k are both proportional to a symmetrical 
effective collision strength, = T^k, which is the source of uncertainty in the collisional 
transition rates. Assuming that the spectral model is arranged in increasing level energy 
order A^ i = whenever k < i. We note that the second term in the denominator of 
the above equation is the inverse of the lifetime of level i, i.e., r, = (5^j<i ^j') -1 - This 
is important because lifetimes are generally dominated by a few strong transitions, which 
are much more accurately determined than the weak transitions. Thus, Tj carries smaller 
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uncertainties than individual A-values. Then, equation (1) can be written as 



Y,k^i N k{ n eqk,i + A k .i) 

n e jn\ + 1 



(2) 



where n| is the so-called critical density of level i and is defined as n\ 



The uncertainty in the population of level i, 5Ni, can be computed as 




(3) 



The first three terms on the right hand side of this equation represent direct propagation of 
uncertainties from atomic rates to or from level i. The last term in the equation correlates 
the uncertainty in level i with the uncertainties in the level populations of all other levels 
that contribute to it . Then, 



where k { = Y<k^i N k( n e<lk,i + A k>i ). 

This linear set of equations yields the uncertainties in the populations of all levels. 
Before proceeding to solve these equations it is worth pointing out some important 
properties: (1) uncertainties are obtained relative to the computed level populations 
regardless of the normalization adopted for these. This is important because while some 
spectral models compute population relative to the ground level other models solve for 
normalized populations such that ^ N k is either 1 or the total ionic abundance. Though, 
the equation above is generally applicable regardless of the normalization adopted. (2) 
In the high density limit, n e oo, the right hand side of the equation goes to zero, 




(4) 
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thus the population uncertainties naturally go to zero as the populations approach the 
Maxwell-Boltzmann values (LTE conditions). (3) By having an analytical expression for the 
propagation of uncertainties one can do a detailed analysis of the spectral model to identify 
the key pieces of atomic data that determine the quality of the model for any plasma 
conditions. (4) The set of linear equations for the uncertainties needs to be solved only once 
for any set of conditions and the system is of the same size as that for the level populations. 
This is unlike Monte Carlo approaches that require solving population balance equations 
hundreds of times, which makes real-time computation of uncertanties impractical. 

The set of equations above can be readily solved by writing them as 

Bx = b, (5) 

where X{ = (5Ni/Ni) 2 , and the matrix and vector elements of B and b are given by the 
equation [2j 

Figure 1 shows the populations and population uncertainties for the first four excited 
levels of O III as a function of the electron density at a temperature of 10 4 K. For 
this computation we have assumed 5% uncertainties in the lifetimes, 10% uncertainties 
in individual A-values, and 20% uncertainties in the effective collision strengths. The 
levels considered here are 2p 2 3 Po,i,2, 1 -D2, and Sq. It is seen that levels 2 through 5 
have maximum uncertainties, ~20%, in the low density limit where the populations are 
determined by collisional excitations from the ground level. As the electron density increases 
thermalization of levels with similar energies and radiative cascades start becoming more 
important, which diminishes the contribution of uncertainties in collision strengths and 
enhances the importance of uncertainties in A-values. For high densities all population 
uncertainties naturally go to zero as the populations approach the Boltzmann limit. 
Another thing to notice is that, the population uncertainties exhibit multiple contributions 
and peaks as the metastable levels 3 Pi and 3 P2 become populated and the uncertainties in 
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these propagate through higher levels. 

Figure 2 shows the populations, relative to the ground level, and population 
uncertainties for the first eight excited levels of Fe II as a function of the electron 
density at a temperatu r e of 10 4 K. For these calculations we use atomic data as in 



Bautista and Pradhan 



(119981 ) and assume uncertainties of 5% in the lifetimes, 10% in 
individual A-values, and 20% in the effective collision strengths. The levels considered here 
are 3ci 6 4s 5 D 9 / 2 j/2, 5/2,3/2 and 3c? 7 4 i 7 9/2,7/2, 5/2,3/2- An interesting characteristic of the Fe 
II system is that the 3c? 7 4 -F9/2 excited level is more populated, at least according to the 
atomic data adopted here, than the ground level at densities around 10 4 cm -3 , typical of H 
II regions. Moreover, under these conditions only ~20% of the total Fe II abundance is in 
the ground level. This means that unlike lighter species, where excitation is dominated by 
the ground level or the ground multiplet, in Fe II all metastable levels are strongly coupled 
and uncertainties in atomic data are expected to propagate in a highly non-linear fashion. 

In Figure 3 we present the population errors for the lowest 52 levels of Fe II at 
T e = 10 4 K and n e = 10 4 cm -3 . These are all even parity metastable levels, except for the 
ground level. The figure shows the total estimated uncertainties together with the direct 
contributions from uncertainties in the collision strengths and A-values (first and second 
terms on the right hand side of equation |2]) and the contribution from level uncertainty 
coupling. It is observed that the collision strengths are the dominant source of uncertainty 
for all levels except level 6 (a 4 F 9 / 2 )- For this level the uncertainty is dominated by the 
A-values and uncertainty couplings with levels of its own multiplet and levels of the ground 
multiplet. This is important because we find that the a A Fg/ 2 level makes the largest 
contribution to the uncertainties in 36 of the lowest 52 levels of Fe II. Unfortunately, the 
atomic data for the 4 -Fg/2 level are among the most uncertain parameters of the whole Fe II 
system, as we discuss in the next section. 
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Fig. 1. — Level populations relative to their Boltzmann limits (upper panel) and relative 
level population uncertainties (lower panel) for the 2p 2 3 Pi (i = 2), 3 P 2 (i = 3), 1 D 2 (i = 4), 
and 1 S'o (i = 5) excited levels of III. 
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Fig. 2. — Level populations relative to the ground level (upper panels) and relative level 
population uncertainties (lower panels) for the 3d 6 As 6 D 7 / 2 (i = 2), 6 D 5 / 2 (i = 3), 6 -D 3 / 2 
(2 = 4), e D 1/2 (2 = 5), 3d 7 4 F 9/2 (2 = 6), 3d 7 4 F 7/2 (2 = 7), 3d 7 4 F 5/2 (2 = 8), and 3d 7 4 F 3/2 
(2 = 9). 
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Fig. 3. — Estimated level population uncertainties for the lowest 52 levels of Fe II at T e = 10 4 
K and n e = 10 4 (black line). Here we assume uncertainties for lifetimes, A- values, and 
collisional rates of 5%, 10%, and 20% respectively. The figure also depicts the contributions 
from uncertainties in collisional rates (blue line), A- values (green line), and coupling of 
uncertainties among all levels (red line). 
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3. Estimating true uncertainties in atomic data 

In the previous section we adopted general uncertainties for lifetimes, A-values for 
forbidden transitions, and effective collision strengths of 5%, 10%, and 20%, respectively. 
In absence of generally accepted procedures to estimate uncertainties in theoretical 
atomic data, these kind of numbers are often cited in the literature as g eneral guidelines ; 



howev er, uncertainty estimates on specific rates are rarely provided. In 



Bautista et al. 



(120091 ) we proposed that uncertainties in gf-values could be estimated from the statistical 
dispersion among the results of multiple calculations with different methods and by 
different authors. The uncertainties can be refined by comparing with experimental or 
spectroscopic data whenever available, although these also have significant associated 
uncertainties. This approach is similar to what has been done for many years by the Atomic 
Spectroscopy Data Center at the Nation Institute of Standards and Technology (NIST; 



http://www.nist.gov/pml/data/asd.cfm) in providing a 'critical compilation' of atomic 



data. 

In estimating uncertainties from the dispersion of multiple results one must keep in 
mind some caveats: (a) Small scatter among rates is obtained when the computations 
converge to a certain value, yet such a convergence is dependant on the maximum size of 
the quantum mechanical representation treatable at the time of the computation. Thus, 
there is no guarantee that every seemingly converged result is indeed correct, as some values 
may result from local minima in the parameter space, (b) Large scatter among different 
calculations is expected in atomic rates where configuration interaction and level mixing 
lead to cancellation effects. The magnitude of these effects depends on the wave-function 
representation adopted. Thus, some computations maybe a lot more accurate than others 
for certain transitions and if we knew which computation is the most accurate, then the 
scatter among all different computations may overestimate the true uncertainty. However, 
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detailed information about configuration and level mixing for every transition is rarely 
available in the literature. Nevertheless, in absense of complete information about every 
transition rate from every calculation, a critical comparison between the results of different 
calculations and other sources of data, if available, provides a reasonable estimate of the 
uncertainty in atomic/molecular rates. 

Are the statistical dispersion values realistic uncertainty estimates? To answer this 
question we look at the intensity ratios between emission lines from the same upper level as 
obtained from observed astronomical spectra and theoretical predictions. The advantage of 
looking at these ratios is that they depend only on the A-values, regardless of the physical 
conditions of the plasma. Thus, the ratios ought to be the same in any spectra of any 
source, provided that the spectra have been corrected for extinction. Fe II yields the richest 
spectrum of all astronomically abundant chemical species. Thus, high resolution optical 
and near-IR [Fe II] lines are the best suited for the present experiment. 

One hundred thirty seven [Fe II] lines are found in the HST/STIS archived spectra of 
the Weigelt blobs of rj Carinae. Six medium dispersion spectra (i?=6000 to 10,000) of the 
blobs were recorded between 1998 and 2004 at various orbital phases of the star's 5.5-year 
cycle. Seventy eight [Fe II] lines are also present in the deep e chelle spectrum (-B=30 000) 



of the Herbig-Haro object (HH 202) in the Orion nebula from 



Mesa-Delgado et al. 



(12009h 



The importance of having multiple spectra from different sources and different instruments 
must not be overlooked. Multiple measurements of the same line ratio minimize the 
likelihood of systematic errors due to unidentified blends, contamination from stellar 
emission, and instrumental effects. 

From the observations, there are 107 line ratios reasonably well measured from the 
spectra. The ratios are defined as 

ratio = max (Fx, F2) jrain (Fx, F2), (6) 
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where F\ and F2 are the measured fluxes of two lines from the same upper level. Here, it 
is important that the minimum of the two fluxes is put in the denominator for the ratio. 
Thus, the line ratios are unconstrained and they are all equally weighted when comparing 
with theoretical expectations. Figure H] illustrates a few line ratio determinations from 
several measurements from spectra of rj Carinae and HH 202, as well as from various 
theoretical determinations. In practice, we perform up to four measurements of every 
observation for different spectral extractions along the CCD and different assumptions 
about the continuum and the noise levels. Thus, we see that the scatter between multiple 
measurements of a given ratio greatly exceed the statistical uncertainties in the line flux 
integrations. Moreover, the scatter between measured line ratios often exceeds the scatter 
between theoretical predictions. Full details about the Fe II spectra and measurement 
procedures will be presented in a forthcoming paper, where we will also present our 
recommended atomic data for Fe II. 

For the present work we consider seven different computations of A-values for Fe II. 



These are the SUPERSTRUCTUR 



Quinet. Le Douneuf. and Zeippen 



and relativistic Hartree-Fock (HFR ) calculations by 



Deb and Hibbert 



(119961 ). the recent CIV3 calculation of 
( 1201 ll ). and various new HFR and AUTOSTRUCTURE calculations that extend over 
previous works. Figure [5] presents a sample of theoretically calculated lifetimes and 
transition yields in Fe II. The yields are defined as yij = A it j x Tj. From the dispersion 
among various results, the average uncertainty in lifetimes for all levels of the 3d 7 and 3d 6 4s 
configuration is 13%. More importantly, it is found that the uncertainty in the critically 
important a 4 -Fg/2 level is ~ 80%, due to cancelation effects in the configuration interaction 
representation of the a 4 -F 9 / 2 — a 6 -D 9 / 2 transition. 

We compared the observed line ratios described above with the predictions from 
different sets of theoretical A-values. Without uncertainty estimates for the theoretical 
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Fig. 4. — Emission line ratios from transitions from the same upper level. The first nine 
points from left to right results from our measured intensities in the HST/STIS spectra of the 
Weigelt blobs of rj Carinae. The tenth point is the measured ratio in the echelle spectrum of 
HH 202. The last point to the right depicts the average of all measurements and uncertainties 
given by the standard deviation. The horizontal lines represent the predictions from several 
different computations of A-values. 
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values, the reduced-^ 2 values from these comparison range from 2.2 to 3100 for the different 
sets of A-values. On other hand, if one adopts average A-values from all calculations 
and uncertainties from the resultant standard deviations the reduced-^ 2 is 1.03. This is 
indicative of well estimated uncertainties, neither underestimated nor overestimated, and 
within these uncertainties there is good agreement between theoretical and experimental 
line ratios. The comparison between observed and theoretical line ratios, including 
uncertainties, is presented in Figure [61 

Figure [7] shows the estimated lifetime uncertainties for the lowest 52 levels of Fe II. 
The figure also presents the level population uncertainties that results from the present 
uncertainties in lifetimes and transition yields for a plasma with T e = 10 4 K and n e = 10 4 
cm" 3 . Here, the adopted uncertainties in the collision strengths are kept at 20% for all 
transitions. By far, the most uncertain lifetime is that of the important a 4 Fg/ 2 level (i = 6), 
yet the way that this uncertainty propagates through level populations depends on the 
density of the plasma. For electron densities much lower than the critical density for the 
level the uncertainty in the lifetime reflects directly on the level population for that level. 
This is seen at n e = 10 4 cm" 3 for levels ~18 and higher. However, as the density increases 
the uncertainties in the level populations become incresingly dominated by the collision 
strengths. This effect is clearly illustrated in Figure |HJ 

4. Uncertainties in Emission Line Emissivities and Diagnostic Line Ratios 

The line emissivity, in units of photons per second, of a transition % — > f, with i > f, is 

j iif = N { x A u . (7) 

In computing the uncertainty in jij one must to account for the fact that iVj and A it f 
are correlated, because the latter appears in the denominator term of Equation (1) that 
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Fig. 5. — Theoretically calculated lifetimes and transiti on yields in Fe II. The calculations 



depict ed are SST: SUPERSTRUCTURE computation by 



jiaad ); h 



Quinet. Le Douneuf. and Zeippen 



results by 



R: HFR calculation also by Quinet et al.; HFRn: our new HFR calculation; CIV3: 



Deb and Hibbert 



(120111 ); ATS21, ATS2, and ATS3: our new AUTOSTRUCTURE 
calculations that extend over Quinet et al. The last point to the right of each panel depicts 
the average value of the various determinations. The uncertainty bars for this point are set 



by the statistical dispersion between all values. 
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Fig. 6. — Line ratios between line from the same upper level measured from optical nebular 
spectra vs. theoretical predictions. 
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Fig. 7. — The upper panel presents the estimated uncertainties in lifetimes for the lowest 
52 levels of Fe II. The lower panel is like Figure [3] but from uncertainties in lifetimes and 
radiative yields estimated from the dispersion among various calculations. 
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Fig. 8. — Like [2] but from uncertainties in lifetimes and radiative yields estimated from the 
dispersion among various calculations. 
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determines JVj. This is important because the most frequently observed lines from any 
upper level are usually those that dominate the total decay rate for the level, i.e., the 
inverse of the level's lifetime. It is convenient to re-write the above equation as 

■ _ A U 

»< Ej %j + Ej A i,i 

Combining this equation with Equation |2] one finds 




This equation can be readily evaluated from the level populations and uncertainties already 
known. The equation has various interesting properties: (1) the equation is independent of 
the physical units used for the emissivities; (2) in the high density limit, as the uncertainty 
in the level population goes to zero, the uncertainty in the emissivity is the same as in the 
A- value. 

Figure M depicts uncertainties in emissivity for a sample of strong IR, near-IR, and 
optical [Fe II] lines. These are computed at 10 4 K. The uncertainties in the collision 
strengths are 20% and the uncertainties in the lifetimes and A-values are those estimated 
in the previous section. The behavior of these uncertainties for different physical conditions 
is complex. Let us look, for instance, at the uncertainty of emissivity of the 5.3/im line 
(a 4 Fg/2 — a 6 Dg/2] 6—^1) whose behavior is contrary to the uncertainty in the population 
of the a 4 -Fg/2 level (see Figured]). According to equations [2] and HI in the low density limit 

3i,f ->• Yl N k n eQk,i ( ^ AlJ A )■ ( 10 ) 

In the case of the a 4 -Fg/2 level the 5.3/xm transition dominates the total decay rate of level 
and the ratio ^j,// Ej<iAi ^ s essentially 1. Thus, the uncertainty in the A^i rate cancels 
out at low electron densities and the uncertainty in the emissivity is small despite a large 
uncertainty in the level population. By contrast, at high densities the population of the 




Fig. 9. — Uncertainties in [Fe II] line emissivities at 10 4 vs. n e . The transitions shown are: (a) 
25.9 //m (a 6 D 7/2 -a 6 D 9/2 ); (b) 5.33 /im (a 4 F 9/2 -a 6 D 9/2 ); (c) 1.256 fj,m (a 4 D 7/2 -a G D 9/2 ); 
(d) 8616.8 A(a 4 P 5 /2-a ^9/2); (e) 7155.2 A(a 2 G 9/2 -a 4 F 9/2 ); (f) 5527.4 A(a 2 £> 5 /2- a 4 F 7/2 ); 
(g) 4889.7 A(6 4 P 5 /2 - a ^7/2); (h) 5261.6 A(a 4 # 11/2 - a 4 F 7/2 ). 
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level approaches the Boltzmann limit and the uncertainty in the emissivity is solely given 
by that in A 6 i, which is ^ 



0. 



A line emissivity ratio between two lines is given by 
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jg,h 
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g,h 



11] 



where AEij is the energy difference between levels % and / and we have used emissivties 
in units of energy per second. In computing the uncertainty in this line ratio one must 
account for the fact that the emissivities are correlated. Moreover, a general expression for 
the uncertainty must account for cases where % — g, in which case the uncertainty in the 
ratio would depend only on the A-values. The uncertainty is the ratio is given by 



if J 
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g,h 



where 



Thus, 
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From Equation |2] we find {dN,/dA gth ) = NiNg/n, for h = i, {dNi/dA g:h ) = -Nf/m for 
g = i, and (dNi/dA g h ) = otherwise. 



In the general case of a ratio involving several lines in the numerator and/or 

kf 



denominator, i.e., 
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(14) 
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Fig. 10. — [Fe II] emissivity line ratios (upper panel) and uncertainties (lower panel) at 10 4 
K vs. electron density. 
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the uncertainty is 




i,f}Hi,f} H{g,h}j{9,h} 
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Figure [TO] shows a sample of line ratios between IR and optical lines and their 
uncertainties. The uncertainties exhibit complex behaviour with changes in density and 
temperatures. In general, line ratios are only useful as diagnostics when the observed ratio 
lies around middle range of the theoretical ratio. Moreover, it is very important to know the 
uncertainties in the ratios when selecting appropriate diagnostics from a given spectrum. 



We presented a method to compute uncertainties in spectral models from uncertainties 
in atomic/molecular data. Our method is very efficient and allows us to compute 
uncertainties in all level populations by solving a single algebraic equation. Specifically, 
we treat the case of non-LTE models where electron impact excitation is balanced by 
spontaneous radiative decay. However, the method can be extended to ionization balance 
and additional excitation mechanisms. 

Our method is tested in O III and Fe II models, first by assuming commonly assumed 
uncertainties and then by adopting uncertainties in lifetimes and A-values given by the 
dispersion between the results of multiple independent computations. Moreover, we show 
that uncertainties taken this way are in practice very good estimates. 

Then we derive analytic expresions for the uncertainties in line emissivities and line 
ratios. These equations take into account the correlations between level populations and 



5. 



Conclusions 
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line emissivities. Interestingly, the behaviour of uncertainties in level populations and 
uncertainties in emissivities for transitions from the same upper levels are often different 
and even opposite. This is the case, in particular, for lines that result from transitions that 
dominate the total dacay rate of the upper level. Then, the uncertainties in A-values for 
the transitions that yield the lines cancel out with the uncertainties in the lifetimes of the 
levels. In terms of emission line ratios, it is also found that knowledge of the uncertainties in 
the ratios is essential selecting appropriate ratios for density and temperature diagnostics. 

At present, we are in the process of estimating uncertainties in atomic data for species 
of astronomical interest. Our uncertainty estimates and analysis of the uncertainties in 
various spectral models, ionic abundance determinations, and dianostic line ratios will be 
presented in future publicaitons. 
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